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Abstract 



o 

■ Microscopic origin of the magnetic field observed in compact stars is studied in 

quark matter. Spontaneous spin polarization appears in high-density region due to 
the Fock exchange term. On the other hand, quark matter becomes unstable to form 
spin density wave in the moderate density region, where restoration of chiral symmetry 
plays an important role. Coexistence of the magnetism and color superconductivity is 
briefly discussed. 



1 Introduction 



Nowadays it is widely accepted that there should be realized various phases of QCD in 
temperature (T) - density (p#) plane. When we emphasize the low T and high pB region, 
the subjects are sometimes called high-density QCD. The main purposes in this field should 
be to figure out the new phases and their properties, and to extract their symmetry breaking 
pattern and low-energy excitation modes there on the basis of QCD. On the other hand, these 
studies have phenomenological implications on relativistic heavy-ion collisions and compact 
stars like neutron stars or quark stars. 

In this lecture we'd like to address magnetic properties of quark matter. We shall see 
various types of magnetic ordering may be expected in quark matter at finite density or 
temperature. They arise due to the quark particle-hole (p — h) correlations in the pseudo- 
scalar or axial-vector channel. We first discuss the ferromagnetic phase transition and then 
another magnetic feature in the second part. 

Phenomenologically the concept of magnetism should be directly related to the origin 
of strong magnetic field observed in compact stars [T]; e.g., it amounts to O(10 12 G) at the 
surface of radio pulsars. Recently a new class of pulsars called magnetars has been discovered 
with super strong magnetic field, B s ~ 10 14 ~ 15 G, estimated from the P — P curve 0. First 
observations are indirect evidences for super strong magnetic field, but discoveries of some 
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absorption lines stemming from the cyclotron frequency of protons have been currently 
reported 

The origin of the strong magnetic field has been a long standing problem since the first 
pulsar was discovered pQ . A naive working hypothesis is the conservation of the magnetic flux 
and its squeezing during the evolution from a main-sequence progenitor star to a compact 
star, B oc R~ 2 with R being the radius. Taking the sun as a typical one, we have B ~ 10 3 G 
and R ~ 10 10 cm. If it is squeezed to a typical radius of usual neutron stars, R ~ 10km, the 
conservation of the magnetic flux gives 10 11 G, which is consistent with the observations for 
radio pulsars. However, we find R ~ 100m to explain B ~ 10 15 G observed for magnetars, 
which leads to a contradiction since the Schwatzschild radius is O(lkm) much larger than 
R. 

Since there should be developed hadronic matter inside compact stars, it would be reason- 
able to consider a microscopic origin of such strong magnetic field: ferromagnetism or spin 
polarization is one of the candidates to explain it. Makishima also suggested the hadronic 
origin of the magnetic field observed in binary X-ray pulsars j^. 

When we consider the magnetic-interaction energy by a simple formula, E mag = \iiB 
with the magnetic moment, /ij = ei/(2m i ) ) we can easily estimate it for B = O(10 15 G); 
it amounts to several MeV for electrons, while several keV for nucleons and 10 - 100 keV 
for quarks. This simple consideration may imply that strong interaction gives a feasible 
origin for the strong magnetic field, since its typical energy scale is MeV. The possibility of 
ferromagnetism in nuclear matter has been elaborately studied since the first pulsars were 
observed, but negative results have been reported so far [3]. In the first part of this lecture, 
we consider its possibility in quark matter from a different point of view (Hj- 

In the second part we discuss another magnetic aspect in quark matter at moderate 
densities, where the QCD interaction is still strong and some non-perturbative effects still 
remain. One of the most important features observed there is restoration of chiral symmetry; 
the blocking of qq excitations due to the existence of the Fermi sea gives rise to restoration of 
chiral symmetry at a certain density and many people believe that deconfinement transition 
also occurs at almost the same time. There have been proposed various types of the p-h 
condensations at moderate densities [7J [H] , in which the p-h pair in scalar or tensor channel 
has the finite total momentum indicating standing waves (the chiral density waves). The 
instability for the density wave in quark matter was first discussed by Deryagin et al. [7j at 
asymptotically high densities where the interaction is very weak, and they concluded that 
the density-wave instability prevails over the BCS one in the large N c (the number of colors) 
limit due to the dynamical suppression of colored BCS pairings. 

In general, density waves are favored in 1-D (one spatial dimension) systems and have 
the wave number Q = 2kp according to the Peierls instability e.g., charge density waves 
(CDW) in quasi- 1-D metals. The essence of its mechanism is the nesting of Fermi surfaces 
and the level repulsion (crossing) of single particle spectra due to the interaction for the 
finite wave number. Thus the low dimensionality has a essential role to produce the density- 
wave states. In the higher dimensional systems, however, the transitions occur provided the 
interaction of a corresponding (p-h) channel is strong enough. For the 3-D electron gas, 
it was shown by Overhauser JU] that paramagnetic state is unstable with respect to the 
formation of the static spin density wave (SDW), in which spectra of up- and down-spin 
states are deformed to bring about the level crossing due to the spin exchange interactions, 
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while the wave number does not precisely coincide with 2kp because of the incomplete nesting 
in higher dimension. 

We shall see a kind of spin density wave develops there, in analogy with SDW mentioned 
above.. It occurs along with the chiral condensation and is represented by a dual standing 
wave in scalar and pseudo-scalar condensates (we have called it 'dual chiral- density wave', 
DCDW). DCDW has different features in comparison with the previously discussed chiral 
density waves [7||Hj. One outstanding feature concerns its magnetic aspect; DCDW induces 
spin density wave. 

2 Ferromagnetism in QCD 
2.1 A heuristic argument 

Quark matter bears some resemblance to electron gas interacting with the Coulomb potential; 
the gluon exchange interaction in QCD has some resemblance to the Coulomb interaction in 
QED , and color neutrality of quark matter corresponds to charge neutrality of electron gas 
under the background of positively charged ions. 

It was Bloch who first suggested a mechanism leading to ferromagnetism of itinerant 
electrons [TT]. The mechanism is very simple but largely reflects the Fermion nature of 
electrons. Since there works no direct interaction between electrons as a whole, the Fock 
exchange interaction gives a leading contribution. Then it is immediately conceivable that 
a most attractive channel is the parallel spin pair, whereas the anti-parallel pair gives null 
contribution (see Eq. (|SJ) below). This is nothing but a consequence of the Pauli exclusion 
principle: electrons with the same spin polarization cannot closely approach to each other, 
which effectively avoid the Coulomb repulsion. On the other hand a polarized state should 
have a larger kinetic energy by rearranging the two Fermi spheres. Thus there is a trade-off 
between the kinetic and interaction energies, which leads to a spontaneous spin polarization 
(SSP) or FM at some density. One of the essential points we learned here is that we need 
no spin-dependent interaction in the original Lagrangian to see SSP. 

Then it might be natural to ask how about in QCD. We list here some features of QCD 
related to this subject. (1) the quark-gluon interaction in QCD is rather simple, compared 
with the nuclear force; it is a gauge interaction like in QED. (2) quark matter should be 
a color neutral system and only the Fock exchange interaction is also relevant like in the 
electron system. (3) there is an additional flavor degree of freedom in quark matter; gluon 
exchange never change flavor but it comes in through the generalized Pauli principle. (4) 
quarks should be treated relativistically, different from the electron system. 

The last feature requires a new definition and formulation of SSP or FM in relativistic 
systems since "spin" is no more a good quantum number in relativistic theories; spin couples 
with momentum and its direction changes during the motion. It is well known that the 
Pauli-Lubanski vector is the four vector to represent the spin degree of freedom in a 
covariant form; the spinor of the free Dirac equation is the eigenstate of the operator, 
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where a 4-axial- vector s.t. 



c + 



k(C-k) 

ra(E k + m) ' 



k-c 



m 



(2) 



with the axial vector £. We can see that is reduced to a three vector (0, £) in the rest 
frame, where we can allocate £ = (0,0, ±1) to spin "up" and "down" states. Thus we can 
still use C to specify the two intrinsic polarized states even in the general Lorentz frame. 
The Fock exchange interaction, /kc,q£'> between two quarks is then given by 



/kf,qC - 9' 



ra m 



E k E q 9m 2 
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in the lowest order, where the spin dependent term renders 

1 



a ■ b = 



rat 



-(k-CKq-O+^C-C 



+ {ra{E k + m) (C • q) (C' • q) + ra(E q + ra) (£' • k) (C ■ k) 
+ (k • q)(C • k)(C • q)} /(E k + ra)(E q + ra)] . 



(4) 



It exhibits a complicated spin-dependent structure arising from the Dirac four spinor, while 
it is reduced to a simple form, 



~9 9 



2 a i + C-C 



q) 



(5) 



in the non-relativistic limit as in the electron system. Eq. (jSJ) clearly shows why parallel spin 
pairs are favored, while we cannot see it clearly in the relativistic expression (0J). We have 
explicitly demonstrated that the ferromagnetic phase should be realized at relatively low 
density region jHj: quark matter is spontaneously polarized at low densities like the electron 
gas. The phase transition is of first order and the critical density is around the nuclear 
density for quark mass ra q ~ 300MeV and coupling constant a c ~ 2.2, which are used in 
the MIT bag model. It should be interesting to see a recent reference j!4j . where the author 
also found the ferromagnetic phase transition at low densities within the perturbative QCD 
calculation beyond the lowest-order diagram. 



2.2 Self-consistent calculation 

If we understand FM or magnetic properties of quark matter more deeply, we must proceeds 
to a self-consistent approach, like the Hartree-Fock theory, beyond the previous perturbative 
argument. 

We begin with an OGE action: 



l int 



(6) 



where D^ u denotes the gluon propagator. By way of the mean-field approximation, we have 

d 4 p 



IMF 



(27T) 



-Mv)G- A \p)^{p). 



(7) 
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The inverse quark Green function G involves various self-energy (mean- field) terms, of 
which we only keep the color singlet particle- hole mean-field V(p), 

Ga(p) _1 =$-m + {i + V{p). (8) 

Taking into account the lowest diagram, we can then write down the self-consistent equations 
for the mean-field, V: 

- V(k) = Hs) 2 J J^u{-iD^(k - P)}l^{-iG A {p)} lv ^. (9) 



V = 



r B 



Nc-l 

Nc 



Figure 1: Graphical interpretations of the equation Q with coefficient in front of R.H.S. 
given by N c . 

Applying the Fierz transformation for the OGE action (jUJ) we can see that there ap- 
pear the color-singlet scalar, pseudo-scalar, vector and axial-vector self-energies by the 
Fock exchange interaction. In general we must take into account these self-energies in V, 
V = U s + jzUps + 7 /i £7£ t + IpH^av w hh the mean-fields C/j. Here we retain only U s , U®, U% v 
in V and suppose that others to be vanished. 

Problem: Show that there is no tensor mean-field as a consequence of chiral sym- 
metry in 

According to the above assumptions and considerations the mean-field V in Eq.fjHJ) ren- 
ders 

V = w s U A , U A = U 3 av , (10) 

with the axial-vector mean-field U a 1 • 

The poles of Ga(p), detG A ^(po=e n )=0, give the single-particle energy spectrum: 

6n = ±e± (11) 

e± = ^P 2 + U 2 A + m 2 ± 2y/m*U A + (p • U A ) 2 , (12) 

where the subscript ± in the energy spectrum represents spin degrees of freedom, and the 
dissolution of the degeneracy corresponds to the exchange splitting of different "spin" states 
[TJ. 

There are two Fermi seas with different volumes for a given quark number due to the ex- 
change splitting in the energy spectrum. The appearance of the rotation symmetry breaking 

1 Since the scalar and vector mean- fields only renormalize the mass and the quark- number chemical 
potential, respectively, we discard them here for simplicity. 
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term, oc p ■ in the energy spectrum implies deformation of the Fermi sea: thus rotation 
symmetry is violated in the momentum space as well as the coordinate space, 0(3) — > 0(2). 
Accordingly the Fermi sea of majority quarks exhibits a "prolate" shape (F~), while that of 
minority quarks an "oblate" shape (F + ) as seen Fig. 2. 






Figure 2: Modification of the Fermi sea as Ua is increased from left to right. The larger 
Fermi sea (F~) takes a prolate shape, while the smaller one (F + ) an oblate shape for a given 
Ua- In the large Ua limit (completely polarized case), F + disappears as in the right panel. 

Here we demonstrate some numerical results; we replaced the original OGE by the "con- 
tact" interaction, D^ v — > —g^/A 2 . Then the self-consistent equation can be written as 



N 2 



-~9 2 



d 3 p ^ af ( ,,U A + s(3 p 



(2vr 



(13) 



within the "contact" interaction, g 2 = g 2 /A 2 . Note that the expression for Ua, Eq. (|T3*jl. 
is nothing but the simple sum of the expectation value of the spin operator over the Fermi 
seas. 



2.3 Phase diagram in temperature-density plane 

We will present the phase diagram in the three-flavor case under mainly two conditions ^H] : 
the chemical equilibrium condition (CEC) p u = fid = Ps and the charge neutral condition 
without electrons (CNC) p u = Pd = Ps, where quark masses are taken as m u = rrid = 5MeV 
and m s 



150 - 350MeV, i.e., p s 



»l,d + m l 



m. 



2 j for T = 0. In both conditions, since 



the spin polarization caused by the axial-vector mean-field is fully enhanced by the quark 
masses, choice of the current quark mass affects the results, especially, the strange quark 
mass is large, and therefore has essential effect on the spin polarization. To get the phase 
diagram, we use the thermodynamic potential Q within the mean-field approximation, 

r d 3 k 

=U,d,! 

d 3 k 



n 



-n c E 

u=±l 



E 



e s (k, rrii,U A 



upi 



+ 1 



s=± 



/d k 
— — -e s (k, m t ,U A ) + 

=«,£t,3 V Z7r ^ 



u 2 

7X 



(14) 



The second term is the vacuum contribution and can be regularized by ,e.g., the proper-time 
method. We can see that the potential Ua reproduces the self-consistent equation Eq. (JT3J) 
in the three-flavor case. 
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Figure 3: Phase diagram for two cases of strange quark masses, 250, 350MeV. 



| Comparison between CEC and CNC| 
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Figure 4: Critical density as a function of g under two conditions: CEC and CNC. 

Fig. El shows the resulting phase diagram under the two conditions mentioned above; 
CEC and CNC. 

The figure shows the critical temperature (the Curie temperature) for a given baryon 
density, CNC tends to facilitate the system being the spin polarization than CEC. This is 
because CNC holds the larger strange-quark density in comparison with CEC. 

Since the axial-vector mean-field comes from the Fock exchange interaction and causes a 
kind of particle-hole condensation, which is enhanced by the Fermi sea contribution, there 
exists a critical density for a given coupling constant g. We show the critical density by 
varying the effective coupling constant g in Fig. 0] The critical density lowers with the 
larger coupling strength, and this tendency is enhanced in the case of CNC. The result also 
indicates that even for the weak-coupling regime in QCD, the spin polarization may appear 
in sufficiently large densities at sufficiently low temperature. 

2.4 Color magnetic superconductivity 

If FM is realized in quark matter, it might be in the color superconducting (CSC) phase [To! 
IT7] . In this section we briefly discuss a possibility of the coexistence of FM and CSC, which 
we call Color magnetic superconductivity 18J. In passing, it would be worth mentioning 
the corresponding situation in condensed matter physics. Magnetism and superconductivity 
(SC) have been two major concepts in condensed matter physics and their interplay has been 
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repeatedly discussed jT^J. Very recently some materials have been observed to exhibit the 
coexistence phase of FM and SC (201 ■ I* 1 our case we shall see somewhat different features, 
but the similar aspects as well. 

We begin with the OGE action (7) and extend it to include the quark pairings. We 
consider the quark pairings on the same Fermi surfaces by using the eigenspinors <f> s , s = 
±1 with the energy (|12|): the pairing function then reads 

A(p) = £ A s (p)B s {p), B a (p) = 7o<Mp)0l(p). (15) 

The structure of the gap function (|15|) is inferred from a physical consideration of a quark 
pair as in the usual BCS theory: we consider here the quark pair on each Fermi surface with 
opposite momenta, p and — p so that they result in a linear combination of J n = 1~. A s 
is still a matrix in the color-flavor space. Taking into account the property that the most 
attractive channel of the OGE interaction is the color antisymmetric 3 state, the quark pair 
must be in the flavor singlet state. Thus we can choose the form of the gap function as 

(A.) fl .. = e a ^A s (16) 

for the two-flavor case (2SC), where a,f3 denote the color indices and i,j the flavor indices. 
In refs. [TH] we can see that the gap functions A s have a polar angle (9) dependence on the 
Fermi surfaces: both the gap functions have nodes at poles (9 = 0, n) and take the maximal 
values at the vicinity of equator (6 = 7r/2), keeping the relation, A_ > A + . This feature 
is very similar to 3 P pairing in liquid 3 He or nuclear matter [2H 122]; actually we can see 
our pairing function Eq. (j!6|) to exhibit an effective P wave nature by a genuine relativistic 
effect by the Dirac spinors. As a consequence, e can say that FM and CSC barely interfere 
with each other [TH] . 

3 Dual chiral density wave 

3.1 Chiral symmetry restoration and Instability of the directional 
mode 

We consider here another type of density wave described as a dual standing wave in the scalar 
and pseudo-scalar mean- fields [2*3*] . It is well known that chiral symmetry is spontaneously 
broken (SSB) due to the quark-anti-quark condensate in the vacuum and at low densities; 
both the scalar and pseudo-scalar densities always reside on the chiral circle of the finite 
modulus, and any chiral transformation by a chiral angle 9 shifts each value on the circle, 
while the QCD Lagrangian is invariant for constant 9. 

The spatially variant chiral angle 9(r) represents the degree of freedom of the Nambu- 
Goldstone mode in the SSB vacuum. The dual chiral density wave (DCDW) is described by 
such a chiral angle 9(r). When the chiral angle has some space-time dependence, there should 
appear extra terms in the effective potential: one trivial term is the one describing the quark 
and DCDW coupling due to the non-commutability of 9{r) with the kinetic (differential) 
operator in the Dirac operator. Another one is nontrivial and comes from the quantum 
effect: the energy spectrum of the quark is modified in the presence of 9{r) and thereby the 
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vacuum energy has an additional term, oc (V#) 2 in the lowest order. This can be regarded 
as an appearance of the kinetic term for DCDW through the vacuum polarization 

3.2 DCDW in the NJL model 

Taking the Nambu-Jona-Lasinio (NJL) model as a simple but nontrivial example, we ex- 
plicitly demonstrate that quark matter becomes unstable for a formation of DCDW above 
a critical density; the NJL model has been recently used as an effective model of QCD, 
embodying spontaneous breaking of chiral symmetry in terms of quark degree of freedom 
|2H1 2 ■ We shall explicitly see the DCDW state exhibits a ferromagnetic property. 
We start with the NJL Lagrangian with Nf = 2 flavors and N c = 3 colors, 

Cnjl = $(0 - m c )4j + G{(^) 2 + (^TstV) 2 ], (17) 

where m c is the current mass, m c ~ 5MeV. Under the Hartree approximation, we linearize 
Eq. (fT7j) by partially replacing the bilinear quark fields by their expectation values with 
respect to the ground state. 

In the usual treatment to study the restoration of chiral symmetry at finite density, 
authors implicitly discarded the speudo-scalar mean-field, while this is justified only for the 
vacuum of a definite parity. We assume here the following mean-fields, 

(ipip) — A cos(q • r) 
(■^ij 5 T 3 ^) = Asin(q-r), (18) 

and others vanish 3 . Accordingly, we define a new quark field ipw by the Weinberg trans- 
formation, 

ipw = exp[ry 5 r 3 q ■ r/2)ip, (19) 

to separate the degrees of freedom of the amplitude and phase of DCDW in the Lagrangian. 
In terms of the new field the effective Lagrangian renders 

C MF = ijj w [i(/> — M — l/2 l5 r 3 ^ w - GA 2 , (20) 

where we put M = —2GA and g M = (0, q), taking the chiral limit (m c = 0). The form given 
in (120)1 appears to be the same as the usual one, except the axial-vector field generated by 
the wave vector of DCDW; the amplitude of DCDW produces the dynamical quark mass 
in this case. We shall see the wave vector q is related to the magnetization: the phase of 
DCDW induces the magnetization. 

2 We can see that the OGE interaction gives the same form after the Fierz transformation in the zero-range 
limit UHI 

3 It would be interesting to recall that the DCDW configuration is similar to pion condensation in high- 
density nuclear matter within the a model, considered by Dautry and Nyman (DN)(2HH2E1j where a and ir° 
meson condensates take the same form as Eq. 1)180 ■ The same configuration has been also assumed for non- 
uniform chiral phase in hadron matter by the use of the Nambu-Jona-Lasinio model However, DCDW is 
by no means the pion condensation but should be directly considered as particle-hole and particle-antiparticle 
quark condensation in the deconfinement phase. 
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Using Eq. (|2~Uj). we can find a spatially uniform solution, ip w = wv^Qo) exp(ip ■ r), with 
the eigenvalues, 



E% = sjEl + |q|V4 ± ^(P • q) 2 + « ^ = (^ 2 + |p| 2 ) 1/2 (21) 
for positive-energy (valence) quarks with different spin polarizations (c.f. (|T2*|) ). 



3.3 Thermodynamic potential 

The thermodynamic potential is given as 
d 3 p 



total 



" ' (2tt) 3 

^val ^va 



[E- - fie- + (e; 

+ M 2 /AG. 
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d 3 p 
(27p 



e~ + e~: 



+ M 2 /AG 



(22) 



where 6± = 6(fi — Ep), /i is the chemical potential and 7 the degeneracy factor 7 = NfN c . 
The first term f2 va i is the contribution by the valence quarks filled up to the chemical po- 
tential, while the second term f2 vac is the vacuum contribution that is apparently divergent. 
We shall see both contributions are indispensable in our discussion. Once fi to tai is properly 
evaluated, the equations to be solved to determine the optimal values of A and q are 



5Q 



total 



5Q 



total 



5A 



Sq 



0. 



(23) 



Since NJL model is not renormalizable, we need some regularization procedure to get a 
meaningful finite value for the vacuum contribution Q vac , which can be recast in the form, 



Z7 



d A p 



trlnS 1 ' 



w ■ 



(24) 



with use of the propagator Sw = (J>— M — l/2T 3 r y 5 ^)~ 1 . Since the energy spectrum is no more 
rotation symmetric, we cannot apply the usual momentum cut-off regularization (MCOR) 
scheme to regularize f2 vac 4 . Instead, we adopt the proper-time regularization (PTR) scheme 
|29j . Introducing the proper-time variable r, we eventually find 



a, 



7 



57T 



3/2 



dr 



T 



5/2 



30 dlh 

00 271 



-(v/pl+A/2+9/2) 2 



+ e 



-(v^f+M^g/2) 2 



n 



ref ; 



(25) 



which is reduced to the standard formula [2Hj in the limit q — > 0. 

The integral with respect to the proper time r is not well defined as it is, since it is still 
divergent due to the r ~ contribution. Regularization proceeds by replacing the lower 
bound of the integration range by 1/A 2 , which corresponds to the momentum cut-off in the 
MCOR scheme. 



4 We would like to stress that the regularization should be, at least, independent of A and q. Otherwise 
it is inconsistent with Eq. I|23|l . 
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Now we examine a possible instability of quark matter with respect to formation of 
DCDW. In the following we first consider the sign change of the curvature of fitotai a t the 
origin (stiffness parameter), (3. Expanding f2 vac with respect to q up to 0(q 2 ), we find 



7A 2 
16tt 2 

nl c + n™ 9 + o(q 4 ), (26) 



fivac = nl c + ^- 2 J(M 2 /A 2 )q 2 + 0(q 4 ) 



where J(x) is a universal function, J(x) = — xEi(— x), with the exponential integral Ei(— x). 
^rac 9 * s the P ure magnetic contribution and provides a kinetic term (oc (V#) 2 ) for DCDW . It 
originates from a vacuum polarization effect in the presence of DCDW and gives a 'repulsive' 
(positive) contribution, so that the vacuum is stable against formation of DCDW. Note that 
it gives a null contribution in case of M = , irrespective of q, as it should be. 

The coefficient of q 2 term in fi™" 9 , the vacuum stiffness parameter, p vac , has a definite 
physical meaning. Since the pion decay constant f n is mass dependent and is given in terms 
of J(x) within the NJL model [2E], flvac can be written as j3 vac = \f 2 . 

For given fi,M and q we can evaluate the valence contribution fi val using Eq. (|2ip. but 
its general formula is very complicated. However, it may be sufficient to consider the small 
q case for our present purpose. Then the thermodynamic potential can be expressed as 
^vai = e V ai(<7) — Wvai(<z)) where e va i(g) and p Y3 i(q) are the energy density and the quark- 
number density, respectively. Expanding e va i(g) up to the 0(q 2 ) we find 



evai(g) ^ 4i + ^M 2 q 2 



— In 



y/fi' 2 - M' 2 M 



+ 0(g 4 ) (27) 



with e° al = 7/8vr 2 [/iv> 2 - M' 2 (2/i 2 - M 2 ) - M 4 ln(/x + yj\i l - M 2 /M)j for normal quark 
matter without DCDW. Thus the valence contribution can be finally written as 

fi val = n° val -^M 2 q 2 H^/M) + 0(q 4 ) 

= ^ai + ^r + 0(g 4 ) (28) 



up to 0(q 2 ), where H(x) = \n(x + Vx^T) and ft° al = e° al -/i P ° al with p° al = ^(^-M 2 f 2 
for normal quark matter. The valence stiffness parameter then reads 

A*i = ~^M 2 H(fi/M) (29) 

Since the function if (x) is always positive and accordingly (3 V(d < 0, the magnetic term VL™°f 
always gives a negative energy and approaches to zero as M — > (triviality). 

We may easily understand why the valence quarks always favor the formation of DCDW. 
First, consider the energy spectra for massless quarks (see Fig. EJ. 

As is already discussed, our theory becomes trivial in this case and we find two spectra 



E? = yJP± + (\Pz\±q/2) 2 , P± = (p*,P«,0), (30) 
which are essentially equivalent to = |p|. 
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-q/2 







q/2 



Pi 



Figure 5: Energy spectra for pj_ = 0. E ± with M = (thick solid and dashed lines). E ± 
with the definite chirality is also shown for comparison (dotted line). We can see there is 
a degeneracy of at p z = for M = 0, while it is resolved by the mass (thin solid and 
dashed lines). 



Problem: Construct the eigenstates of chirality 75 by taking the linear combination 
of tpw- Their spectra are given by E^ 1 = Jp± + (p z ± q/2) 2 . 

There is a level crossing at p = 0. Once the mass term is taken into account this 
degeneracy is resolved and the exchange splitting appears there. Hence it causes an energy 
gain, if q = 0(2 fx); we can see that this mechanism is very similar to that of SDW by 
Overhauser |HJ [TU] . 

Using Eqs. (J22J), (EHJ), (|28jl we write the thermodynamic potential as 



with the total stiffness parameter (3 = (3 vac + f3 va i and the usual NJL expression without 
DCDW, tt NJL = fi° ac (M) + fi vai( M ) + M 2 /AG. The dynamical quark mass M is given 
by the equation, dVt tota \/dM = 0; At the order of q° the dynamical quark mass M° is 
determined by the equation, dtt NJL /dM\ M0 = 0. Since M — M° = 0(q 2 ), DCDW onsets 
at a certain density where the total stiffness parameter f3 becomes negative: the critical 
chemical potential fi cr is determined by the equation, 



In the limit M° — > 0, we find /2 cr = l/2e _TB//2 ~ 0.375... with 7^ being Euler's constant and 
fl = fi/A in the PTR scheme. Note that this is only a sufficient condition for the existence 
of DCDW phase, and we can never exclude the possibility of the first order phase transition 
or metamagnetism Actually, we shall see that DCDW occurs as a first-order phase 

transition. 

3.4 First-order phase transition 

The magnitudes of M and q are obtained from the minimum of the thermodynamic potential 
(j22j) for T = 0. Fig. IH1 shows the contours of Q tot in M-q plane as the chemical potential 



n tom = n NJL + (3q 2 + 0(q 4 ) 



(31) 




(32) 
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(a) fi=OA\ {b)/i=0.48A {c)^=/4,=D.48931A 




M/A M/A M/A 

Figure 6: Contours of fitotai at T = are shown in M — q plane as the chemical potential 
increases, (a) — > (f). The cross in each figure denotes the absolute minimum. 



increases, where the parameters are set as GA 2 = 6 and A = 850 MeV, which are not far 
from those for the vacuum (/x = 0) 

The crossed points denote the absolute minima. There are two critical chemical potential 
/x = fi c i, /x c2 : for the lower densities (Fig. Eta)- (b)) the absolute minimum resides at the point 
(M 7^ 0, q = 0) indicating the SSB phase. At /x — \i c \ (Fig. EJc)) the potential has the two 
absolute minima at (M ^ 0, q = 0) and (M ^ 0,g ^ 0), showing the first-order transition 
to the DCDW phase which is stable for /i cl < /x < /x c2 (Fig. ((Hj)d-e). At /x = /x c2 (Fig. Eff)) 
the axis of M = and a point (M ^ 0,g ^ 0) become minima, the system undergoes the 
first-order transition again to the chiral-symmetric phase. 

Fig. shows the behaviors of the order-parameters M and q as functions of /x at T = 0, 
where that of M without DCDW is also shown for comparison. It is found from the figure 
that DCDW develops at finite range of /x (/x cl < /x < /i C 2)> where the wave number q increases 
with fi, which value is smaller than twice of the Fermi momentum 2/c^(~ 2/x for free quarks) 
due to the higher dimensional effect; the nesting of Fermi surfaces is incomplete in the 
present 3-D system. Actually, the ratio of the wave number and the Fermi momentum (at 
normal phase q = M = 0) becomes q/kp = 1.17 — 1.47 for the baryon-number densities 
p b /p = 3.62 - 5.30 where the DCDW is stable (see Fig. 8). 



3.5 Magnetic properties 

The mean-value of the spin operator is given by 

_ 1 f l q/2±0 

s z = -u ] w ^zUw = 77 — h vac, (33) 

2 2 bp 

with f3 = ^Jpl + m 2 , where "vac" means the vacuum contribution. First note that the integral 
of s z over the Fermi seas should be proportional to q, and the solution with q ^ seems to 
imply FM. However, we can show that PTR gives the vacuum (the Dirac sea) contribution 
oppositely to cancel the total mean-value of the spin operator, which is consistent with 
Eq. (J23J). Instead we can see that the magnetization spatially oscillates, 

M z = (q<7 12 q) = (70CX12) cos(q • r), (34) 
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< 0.4 - 
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_ M/A(q=0) 
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Figure 7: Wave number q and the dynami- 
cal mass M are plotted as functions of the 
chemical potential at T = 0. Solid (dotted) 
line for M with (without) the density wave, 
and dashed line for q. 




0.50 



0.55 



0.60 



fl/A 



Figure 8: Baryon number density func- 
tion of /i. po — 0.16fm~ 3 : the normal nuclear 
density. 



with 



(700-12) 



d 3 p 



2M 



F+-F- (2tt) 3 l M 2 + p 2 



(35) 



which means a kind of spin density wave [30 . 



3.6 Phase diagram at T — fi plane 

To complete the phase diagram we derive the thermodynamic potential at finite temperature 
in the Matsubara formalism. The partition function for the mean-field Hamiltonian is given 
by 



DipDip exp J dr j d 4 r <j r it) ■ M cxp (/-. -q • r) 



i> - 



4G 



n {^ n +^-EUk)} NfNc 

k,n,s=± 



x exp 



1H 



(36) 



where (3 = 1/T, d = — 70^+27 V and uj n the Matsubara frequency. Thus the thermodynamic 
potential VLp is obtained, 

Qp(q,M) = -TlogZp(q,M)/V 

= —NfN c J 7 0E{ ri °g[ e " /3(Bs(k) " M) + 1 ] [e^ E ^ )+ ^l]+E s (k)}+^.(37) 

From the absolute minima of the thermodynamic potential (|37j) . it is found that the 
order parameters at T 7^ behave similarly to those at T = as a function of //, while 
the chemical-potential range of the DCDW at finite temperature, fi c i{T) < \i < /i C 2(T), 
gets smaller as T increases. We show the resultant phase diagram in Fig. Q, where the 
ordinary chiral-transition line is also given. Comparing phase diagrams with and without q, 
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Figure 9: Phase diagram obtained from the thermodynamic potential (|37jl . The solid 
(dashed) line shows the chiral-restoration line in the presence (absence) of DCDW. The 
shaded area shows the DCDW phase. 

we find that the DCDW phase emerges in the area (shaded area in Fig. (|9~|)) which lies just 
outside the boundary of the ordinary chiral transition. We thus conclude that the DCDW 
is induced by finite-density contributions, and has an effect to expand the chiral-condensed 
phase (M 7^ 0) toward low temperature and high density region. 

4 Summary and Concluding remarks 

We have seen some magnetic aspects of quark matter: ferromagnetism at high densities 
and spin density wave at moderate densities within the zero-range approximation for the 
interaction vertex. These look to follow the similar development about itinerant electrons: 
Bloch mechanism at low densities and spin density wave at high densities by Overhauser. 

With the OGE interaction, we have seen ferromagnetism in quark matter at low densities 
(Sect. 2.1) . It would be worth mentioning that another study with higher-order diagrams 
also supports it. These studies suggest an opposite tendency to the one using the zero-range 
interaction (Sects 2.2,2.3). Note that we can also see the same situation for the itinerant 
electrons; the Hartree-Fock calculation based on the infinite-range Coulomb interaction fa- 
vors ferromagnetism at low density region, while the Stoner model, which introduces the 
zero-range effective interaction instead of the Coulomb interaction, gives ferromagnetism 
at high densities. So we must carefully examine the possibility of ferromagnetim in quark 
matter by taking into account the finite-range effect. 

We have seen that dual chiral desnity wave (DCDW) appears at a certain density and 
develops at moderate densities (Sect. 3). It occurs as a result of the interplay between 
the qq and particle-hole correlations. The phase transition is of weakly first order, and the 
restoration of chiral symmetry is delayed compared with the usual scenario. Note that there 
is no soft pions in this phase, and phason may appear istead as a Nambu-Goldstone boson. 

For the discussion of DCDW, we have seen the remarkable roles of the Fermi sea and the 
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Dirac sea: the former always favors DCDW, while the latter works against it. The similar 
situation also appears about the magnetic property of quark matter. The mean value of 
the spin operator over the Fermi seas of valance quarks always gives a finite value in the 
presence of DCDW, which is a kind of ferromagnetism, but the vacuum contribution given 
by the Dirac seas completely cancels it. As a result there is no net spin polarization in this 
case, but we have seen magnetization spatially oscillates instead (spin density wave). This 
is one of the typical examples in which the nonrelativistic picture is qualitatively different 
from the relativistic one by the vacuum effect. 

It would be ambitious to give a scenario based on ferromagnetism of quark matter, which 
can explain the hierarchy of the magnetic field observed in three classes of neutron stars, 
magnetars, radio pulsars and recycled millisecond pulsars. It would be also interesting to 
study some implication of spin density wave or phason on the magnetic properties of compact 
stars. 
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